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Abstract. We analyse the surface density of very faint galaxies 
at the limit of the sky background noise in the field of the 
cluster of galaxies C10024-I-1654. The radial variation of their 
number density in the magnitude bins B = 26 — 28 and / — 
24 — 26.5 displays an (anti)bias magnification effect similar to 
the one observed in A1689 by Broadhurst (1995) for / < 24. 
The study of this effect provides a determination of the radial 
distribution of critical lines of the gravitational lens from which 
we infer the redshift range of the populations seen in B and I. 

The surface density of B galaxies drops rapidly to a well- 
defined inner critical line near the large triple arc. The deple- 
tion extends from Rb ~ 30" to Rb = 53" and the shape of the 
curve can almost be reproduced with two redshift populations 
selected by the two fihers. With our data 60% ± 10% of the B- 
selected galaxies are between ^ = 0.9 and z — 1.1 while most 
of the remaining 40% galaxies should be at a redshift close to 
z = 3. 

The / selected population observed with the lens has a sim- 
ilar but broader depletion with a minimum extending from the 
B inner critical line to Ri = 60". Whatever the cosmological 
model, the very faint /-selected galaxies spread up to a larger 
redshift with about 20% above z > 4. The fact that many 
faint I selected galaxies are not detected in B also favour the 
existence of a more distant population with a redshift range 
between z = 3 and z — 4.5. 

Using a model for the gravitational potential derived from 
a study of the large triple arc seen around the cluster, the loca- 
tions of the two extreme critical lines for the B and I selected 
galaxies favour f2A-dominated flat universes with a cosmologi- 
cal constant ranging from 0.6 to 0.9. The result is confirmed by 
a preliminary investigation of the Broadhurst's effect in A370. 

However, ultra-deep detection of faint distant galaxies 
down to the noise level are technically very difhcult. In this 
first paper we mainly discuss the method to search the so- 
called "last critical line" and several possible effects which may 
bias the results on cosmological parameters. We conclude that 
the systematic measurement of this outer critical line of the 
faint / selected galaxies population around many clusters with 
gravitational arcs of known (low) redshift may help to count 
the number of faint galaxies at very large redshifts beyond the 



possibility of any spectroscopic survey, and eventually to settle 
the issue of the existence of the cosmological constant. 
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1. Introduction 

Broadhurst et al. (1995) and Broadhurst (1995) pointed out the 
interest of measuring the apparent variation of the surface den- 
sity of background galaxies around massive clusters of galaxies. 
This variation is due to competing gravitational magnification 
and deviation effects. It depends on the slope of the galaxy 
counts as a function of magnitude and on the magnification 
factor of the lens which is responsible for a radial stretching 
of the field of view around the cluster. The effect can be used 
to map the region where the magnification factor is maximum 
(critical lines) and to evaluate the mass of the cluster. More 
generally, it provides a way to break the intrinsic degeneracy of 
the inversion methods in studies of gravitational weak shears 
(Schneider & Seitz 1995) . Combined together, these tools pro- 
vide a unique way to recover the total mass and the mass profile 
of gravitational systems. 

The radial surface density of background galaxies up to a 
magnitude limit m can be expressed as 



2.57-1 



N{< m,r) = 7Vo(< m) /^(r) 
where 7 is the intrinsic count slope 
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dlogN{< m) 
dm 
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/i(r) is the magnification factor of the lens and No{< m) the 
intrinsic counts in the absence of a lens, which is obtained 
from counts in a nearby empty field. It is clear from expres- 
sion (1) that a radial magnification bias 7V(< m,r) shows up 
only when the slope 7 is different from the value 0.4; oth- 
erwise, the increasing number of magnifled sources is exactly 
cancelled by the apparent field dilatation and there is no ef- 
fect on N{< m,r). A radial amplification bias should not be 
observed behind a nearby cluster in B{< 26) since the slope is 
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almost this critical value (Tyson 1988). It is only detected in 
the R or I bands for which the slopes are close to 0.3 (Small et 
al. 1995). 

When 7 < 0.3 a decrease of the number of galaxies is 
expected in regions of magnification. The effect is strongest 
at the critical radius around the cluster corresponding to the 
most probable rodshift of the background sources. Since the 
critical radius increases with redshift, two populations of dif- 
ferent redshift will show two different radial depletion curves, 
with their minima deferring by an amount which depends on 
the redshift difference and the cosmological model. Therefore, 
from the study of the radial depletion, one can get informa- 
tion on the redshifts of the background galaxies. Alternatively, 
Broadhurst (1995) emphasised that this effect may be used to 
constrain the geometry of the Universe, provided we know the 
redshift of the two populations. 

A difficulty comes from the fact that the critical radius 
writes rcr oc Dls{zi, Zs,i1m,^A)/ Dos{zi, Zs,flm,iiA) and de- 
pends both from the redshift and cosmology (Dls and Dql are 
respectively the angular distances from the lens to the sources 
and from the observer to the sources). In this paper we show 
that it is possible to disantangle the redshift distribution from 
the geometry of the universe provided the properties of the lens 
are known from the redshift of one critical line, such as the red- 
shift of a giant arc. In the case of C10024-I-1654 we show that a 
population of faint blue galaxies exists that traces a reference 
critical line at low redshift. 

There are only a few clusters for which the properties of 
the lens are known, around which deep B and I observations 
are available on a large field. An analysis of a set of rich lens- 
ing clusters is in progress (Mellier, Fort & Dantel-Fort 1996, 
hereafter paper II). In this paper, we present the results on 
C10024-I-1654 and some preliminary results on A370 in order 
to focus on the basic principles of the method. 

In section 1 we show that we can actually detect in deep 
CFHT images very faint galaxies at the limit of the noise in 
the respective bins B = 26 — 28 and 7 = 24 — 26.5 and that the 
slope of count rates in both magnitude range are appropriate to 
produce an observable magnification bias. This approach has 
been motivated by the recent works done by van Waerbeke et 
al. (1996a) who used the pixel-to-pixel autocorelation function 
(ACF) of the sky background for measuring the shear field 
around lensing clusters. Though the ACF does not need ob- 
ject detection and cannot provide simply informations on the 
galaxy number counts, from the strong signal they measured in 
C10024-I-1654, we can predict that the number density of faint 
sources is high and could be useful for the magnification bias. 

The second section discusses the depletion around 
C10024-I-1654 (z = 0.39). This cluster is ideal because giant 
multiple arcs have been observed which provide consistent lens 
modelling (Kassiola et al. 1992, hereafter KKF; Bonnet et al. 
1994, Walhngton et al. 1995). 

The last section discusses the results including some results 
obtained on A370. 

2. Deep galaxy counts at B > 26 and J > 24 

Deep observations of the cluster 0024-1-1654 (z = 0.39) and of a 
nearby reference field offset by about one degre from the cluster 
center were secured at the prime focus of the CFHT during the 
two nights in B and I with a constant seeing (0.7" in B, 0.55" 
in 7). The total exposure time are 10800 s for each filter. In 



B the field of view is 4.15'x3.1' for the cluster fields and only 
2.2' X 3.5' for the blank field. In 7, the field of view is 7'x7' for 
the cluster and the blank field. The 7-blank field observation 
was contaminated by a strange meteoritic shower (Jenniskens 
et al. 1993) and the corresponding counts is not considered. 
The I count is instead extracted from the Hubble Deep Field 
study. A description of the observations, with a modelling of 
the large triple arc around the cluster, can be found in KKF. 
The ellipticity of the potential is small and we assume that it 
does not change appreciably within 30" to 60" from the cluster 
center. 

We detected very faint galaxies by using the SExtractor 
software (Berlin and Arnouts 1995). The optimization of the 
input parameters necessary to detect extremely faint objects 
was performed from simulations and are available on request 
with a reference blank field in B. The simulated galaxy counts 
parameters of the observations were chosen to be similar to 
those in our CFHT CCD images. 
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Fig. 1. Cumulative B galaxy counts in the fields C10024+1654 (full 
line), in the CFHT blank field (dashed line) and obtained from the 
HDF (dotted line). The agreement between the CFHT blank field 

and the HDF is excellent in all the magnitude bins. The thin line 
with open circles show the cluster members as expected from the 
subtraction of the cluster field to the blank field. The flat shape 
beyond B = 25.5 show that the contamination is negligible in the 
magnitude range 26 < -B < 28 discussed throughout this paper. The 
dotted-dashed curves show an estimate of the fraction of spurious 
detections inferred from the optimized parameters used in SExtrac- 
tor on simulated CCD images. 

One of the main concerns is the significance of detection at 
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faint magnitudes. Since a wrong estimate of the count slopes 
could change the estimate of the magnification, wc performed 
additional simulations of a blank field, free from any objects 
(stars or galaxies), with only a pure poissonian sky background. 
We then optimised the parameters defining galaxies so that 
only a small fraction (< 10%) of spurious events are detected. 
Indeed the detection depends on the noise amplitude and the 
method shall be used on images with good flat-fielding and 
stable noise statistic. Any area with an enhancement of the 
background like scattering of light around bright stars or bright 
galaxies will change the limit of detection locally and shall 
be masked on the image. Finally, the surface number density 
inferred from the galaxy counts in C10024+1654 fields has been 
corrected from the effect of sky occultation by bright cluster 
members. 

Figure 1 shows the B counts in the cluster field and the 
blank field. The counts have been done independently in the I 
filter also. The detection limits are respectively 28 in B, and 
26.5 in /. The number of expected spurious detections coming 
from simulations is also indicated. A plot of the recent Hubble 
Deep Field counts (Couch 96, private communication) is also 
given. The agreement between the HST counts and the CFHT 
counts for the blank field is excellent, giving confidence in our 
method. This is a good exemple of what can be done from the 
ground with good seeing. 

We find slopes of 0.17 ± 0.02 and 0.25 ± 0.03 respectively 
for the counts in B at 26 < B < 27.5 and in / at 24 < / < 26.5 
in the blank fields. These values are close from the HDF counts 
(Table 1) and the Tyson's counts (1988) in B, and compatible 
with the deep I band counts of Small et al (1995). One can see 
on Figure 1 that at B > 26 the contamination by faint cluster 
members is negligible. 

3. The magnification bias in C10024+1654 and the red- 
shifts of the blue and red background sources 

Since the galaxy counts in both B and I have slopes below 0.25 
we can expect a significant depletion of galaxies at some radii 
around the cluster field. The radial surface density obtained 
for B and I are given on figure 2 and 3 respectively. Both 
curves decrease with radius down to a minimum, and increase 
towards A^o at large radius. The B shows a steeper descent. Its 
minimum is pronounced at about 30 arcsec. and extends to 50 
arcsec. from the cluster center. The I curve shows a shallower 
broad dip ranging from 30" to about 60" . As the B and I count 
slopes arc almost identical, the width of the depletion probe the 
relative redshift distribution of the two samples. The difference 
observed beyond 50 arcsec. already shows qualitatively that a 
fraction of the / galax;ies must be at a redshift significantly 
higher than the B sample. IncidentUy we note that the behavior 
of the B-selected galaxies not detected in I is similar to the 
whole sample of B galaxies. 

The depletions of the curves given in figures 2 and 3 can be 
reproduced from a least square fit of the depletion generated 
by a population of galaxies located at various redshift: 

M 

N{< m,r) = Y, Noi{< m) Mi(r-)'-'^*-' . (3) 

where Noi, Hi and 7, are respectively the number density of 
galaxies, the magnification and the slope of the galaxy counts 
of the population located at redshift Zi. The free parameters are 
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Fig. 2. Depletion curve obtained in C10024+1654 from galaxy 
counts in the range 26 < -B < 28 for the B-selelected galaxies. The 
data are the filled circles with error bars. We see that a depletion is 
detected, with a minimum starting at Rb = 30 arcsec. The full line 
show a fit with the model and core radius obtained by KKF. 



the slope of the galaxy counts, which at the starting point of the 
minimisation are assumed to be the slopes of the blank field, 
the redshift and the number density of the sources. Though 
we allow the redshift distribution of galaxies to be described 
with a maximum M of 5 redshift bins (equivalent to 5 critical 
lines), we searched for solutions using a minimum of lens screen 
to describe the distribution of sources in redshift which are 
selected through our two filters. The coefficients Aoi in Table 
1 correspond to the galaxy number density at redshift Zi. It is 
important to note that since the apparent depletion of galaxies 
only lies in the innermost region of the Icnsing cluster, the 
choice of the analytical function representing the mass profile 
is not critical for radius in the range 25 to 150 arcseconds 
(Bonnet et al. 1994). Therefore, we fit the data by using the 
isothermal model of KKF with a core radius Vc = 15 arcsec 
and a velocity dispersion of u = 1300 km/ sec. 

In the following the position of the two extreme critical 
lines are the important output of the fit rather than the rela- 
tive redshift abundances of distant galaxies, a topic that would 
need more additional data on many clusters. Hence we have 
mostly explored the possible range of variation for the radius 
of the extreme critical lines in order to estimate a possible 
dependance on the cosmology. 

For the B galaxies, the best fit of the depletion curve can 
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Table 1. Summary of the galaxy counts and the magnification bias observed in the lensing cluster C10024+1654 in the B and I 
bands. Notot is the total galaxy number density in the cluster field. Nhb, Noh arc the true background components in the two 
bands, and the foreground or cluster members superimposed to the lensed population. The densities are given in gal /arcmin^ 
and are corrected from the obscurations by bright foreground galaxies. The two last columns give the slopes of the galaxy counts 
as inferred from the HDF and from the best fit of the magnification bias. The number density inferred when the depletion is 
reproduced from a best fit is given for the various redshift bins found for the distribution of galaxies. The fit give Ri, Nou, and 
the slope 7^. The rodshifts of each bin docs not depend much on the cosmology and according to the discussion in paragraph 4 
are given for a flat universe Qm = 0.35 and Qa = 0.65. 



be modelled by using only three (almost 2) redshift bins for 
the B sources population, each of them corresponding to a 
slope 7 — 0.20 ± 0.02 in agreement with the slope derived 
from the blank field. The results of the fit is given in table 1 
(reduced = 0.66). The inner critical radius is found at Rb = 
29.8" ±0.4" . This is a robust value which is well constrained by 
the steep shape of the depletion as the radius increases from 
the cluster center. If we allow the core radius to vary around 
the KKF value we find a very similar but marginally best fit 
(X^ = 0.63) with a core radius of 13" and Rb = 30.2". The 
uncertainties coming from a small residual error on the value 
of the core radius of the KKF model is negligible and in the 
following we keep the KKF value. But it is important to note 
that both the core radius and ellipticity of such an isothermal 
model could be determined from a mapping of the anti-bias 
magnification effect {N{x, y; < m)) over the whole CCD image 
(see paper II). 

Whatever the cosmology the KKF model imposes that the 
nearest B galaxies have a redshift between 0.9 and 1.1. In fact 
this value was not unexpected. First, Mellier et al. (1991) ob- 
tained spectra of the images of the arcs whose postions are 
close to Rb = 35" (see figure 2). They are blue objects and 
they should belong to the B population. The presence of bright 
starburst spots clearly visible in the HST images of these arcs 
(Colley et al. 1995) is a strong argument that star formation is 
ongoing in this galaxy. Although a definite redshift cannot be 
given, the spectra provide compelling evidence that it could be 
close to one, in any case larger than 0.88: otherwise we should 
see a proheminent [OII]A3727 emission line in the observed 
spectral range covered by the data. Second, with the shear 
profile in C10024-I-1654, Bonnet et al. (1994) infer that the B 
galaxies responsible for the shear pattern should have an aver- 
aged redshift between 0.9 and 1.2. Third, the redshift survey of 
giant arcs (see Fort and Mellier 1994, Soucail 1995) and arclets 
(Bezecourt and Soucail 1996, Pello et al. 1996), and the lensing 
inversion in A370 and A2218 as well (Kneib et al 1994, 1996) 



provide a fairly large sample of blue-selected galaxies with B 
magnitude ranging from 24.5 to 27, which all have a median 
redshift of 0.85. Most of them are between 0.65 and 1.10, but 
the lower redshifts of arclets have been found behind the lower 
redshift clusters (zi < 0.25) as expected from the properties of 
gravitational lensing. Besides, the average magnitude of these 
samples is about one magnitude brighter (26-27.5 instead of 
24.5-27) so we can reasonnably expect that our faint B galax- 
ies should be somewhat more distant than z = 0.75, whatever 
the evolution model. 

The depletion curve of the B galaxies raises again at 
R3 = 53" corresponding to a redshift of about z = 3.0 ± 0.5 
well determined from the best fit. The redshift error is mainly 
dominated by the uncertainty on the redshift of the inner crit- 
ical line. An important point is that the position at which the 
B depletion curve starts to increase does not depend on the 
number of source screen (Table 1), but only depends on the 
largest redshift screen. Therefore, the depletion curves has the 
remarkable property to provide the spread of the redshift dis- 
tribution, from the width of the depletion, and the highest 
redshift of the sample from the comparison of the radius of the 
upper critical line with a lower critical line. At this stage and 
before the detailed discussion of extreme redshift in paragraph 
4 let's assume that the blue population range from zb = l.±0.2 
to Zb = 3, with a relatively well defined lower limit from the 
convergence of the lens. 

For the faint / galaxies the observed rmmbcr density 
Ni {25 < I < 26.5) is depleted in a larger range of radii. The 
beginning of the dip for the I galaxies is almost the same as 
for the B galaxies, thus their averaged redshift is the same 
(Table 1). This is not surprising as this is the lowest redshift 
for which the gravitational magnification is substantial. The 
minimum of the depletion is broader than for the B-selected 
sample but it has a bump between the two extreme radii. In 
figure 3, the best fit of the data uses the same lens model as for 
the B galaxies, but needs a model distribution of galaxies with 
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4. A tentative to constrain the cosmological parame- 
ters with CI0024+1654 

4-1. The relative scaling method 

The beginning of the dip of the B depletion and the radial 
position at which the number density of I galaxies raises again 
define two extreme critical lines, Rmin and Rmax- Their ratio 
writes: 



Rmaa 



1/2 



(4) 
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Fig. 3. Same plot as figure 2 for the / sample. The width of the 
depletion spread on a larger area. The minimum radius is almost 
similar to the B galaxies. The Jast critical radius is defined as the 
position where the depletion curve starts to increase. This corre- 
sponds to Ri = 60" aJTcsec. 



4 discrete different z corresponding to critical lines Ri. The 
reduced is very good (0.57) if we substract to the data a 
uniformly distributed component corresponding to 14% of the 
count in the blank field. The origin of this small count excess 
and bump in the bottom of the depletion is not clearly under- 
stood so far (cf paragraph 5). If the bump and the relatively 
deep dip near the last critical line are real (Figure 3) they arc 
not perfectly reproduced with the model (table 1). Instead the 
radius of the two extreme critical radii are robust and the best 
fit gives: R2 = 32" ± 1" and R5 = 60" ± 2". Prom the mod- 
elling of the flat depletion curve we conclude that the redshifts 
of the / galaxies are 1. < z/ < 4.5, but with a significant frac- 
tion at large rcdsliift (table 1). According to the KKF model 
the most distant critical line corresponds to an abundance of 
galaxies of 20% at z > 4. This is an extreme critical line where 
the effects of the cosmological parameters can be seen. In fact 
C10024-I-1654 is an ideal configuration where the role of the 
cosmology can bo observed directly from the radius of an ex- 
treme critical line, though its effect is mixted with the redshift 
distribution of the sources. In the next section we consider in 
detail the possibility to obtain simultaneously the redshift dis- 
tributions and the cosmological parameters from the positions 
of critical lines. 



and provides an interesting observable parameter based on the 
two well-measured positions of the smaller and the larger crit- 
ical radii. The dependence on a is weak and is non-existent for 
a singular isothermal sphere. The dependence with Vc is gen- 
erally weak for simple potential since the cluster core radii are 
small with respect to the position of the external critical line 
(see Fort & Mellier for a review). However, it depends on zb, 
zi and the geometry of the universe defined by flm and Oa. 

It is possible to map the volume space (z min , Zmax , , 
in order to find the solutions which gives a Rmax /Rmin in 
agreement with the observations. The results are shown on fig- 
ure 4 for flat universes, and 011 figure 5 for universes with a cos- 
mological constant equal to zero. We have kept all the accept- 
able solutions, with extreme error bars {RBmin = 30" ± 0.5", 
= 60" ±2"), which corresponds to 1.9 < R max I Ftmin < 
2.1, for an isothermal models with a core radius corresponding 
to the best fit parameters found by KKF. The plots only show 
models with zb = 0.85, 0.95 and 1.1, because they are suffi- 
cient for the conclusions. Both figures show that we cam find ac- 
ceptable solutions for amy cosmological models. However, there 
are clear trends for the redshifts of the sources. Whatever the 
acceptable models, they all imply that the clostest / and the B 
galaxies are at redshift lower than 1.0, and the most distant / 
galaxies are at redshift higher than 2.5. A limit of the / galax- 
ies {zi < 6.5) is given by the filter itself which imposes that 
these galaxies have the Lj/a break below 9000A. A lower limit 
at z — 0.85 is given by the spectroscopic informations from the 
giant arc. 

The relative scaling method gives informations on the red- 
shift distribution of galaxies that are almost independant from 
the lens modelling. But it does not constrain much the cosmol- 
ogy as far as the information on the exact radial location of 
the critical lines is not used. These positions depends on the 
lens model through the parameter a. If cr is determined from 
the redshifts of arcs located at low redshifts, the location of 
the most distant / critical line becomes a relevant observation 
for the cosmology. 

4-2. Scaling the critical line with the lens modelling 

From the last critical line Ri which lies at a large distance from 
the cluster center it is possible to compute the angular distance 
ratio Dls/Dos with the lens model. For simplicity, let us just 
express it for the case of an isothermal sphere with a core 
radius rc- The relation between Dls/Dqs and the position of 
the critical line Ri is: 
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Fig. 4. Ratio Ri/Rb versus redshift of the / sources for three red- 
shifts of the blue populations {zb) corresponding to the beginning 
of the dip in the B depletion curve, and various flat universes. The 
horizontal lines give the same quantity as inferred from the obser- 
vations of Rb = 30 ± 0.5 and Ri = 60 ± 2 and computed with a 
singular isothermal sphere with a core radius of 15 arcsec and a ve- 
locity dispersion of 1300 km/scc. The central line is the central value 
and the two other lines show the position at ilcr. Since we only use 
ratios of angular distances, it does not strongly depends on the lens 
modelling. We see that the data are compatible with any models, 
but the redshift of the blue population must be below zb = 0.95. 
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Fig. 5. Same plot as figure 4 but for cosmological models with a 

cosmological constant equal to zero. Again any value of O m is com- 
patible with the observation, provided the redshift of the blue pop- 
ulation observed at the beginning of the dip is lower than 0.95. 



For a given tr, this ratio becomes significantly dependent on 
the cosmological model when the source redshift is large, as 
it is the case for the / population. Hence we search in which 
redshift range and for which set of cosmological parameters, 
the ratio DLs[zi,VLm, Qa)/Dos{zi, ^m, ^a) is equal to the ob- 
served value. 

Figure 8 shows the relation Dls/Dqs to the redshift of 
the sources, for a lens at zi = 0.39. The curves correspond 
to various cosmological models with a zero cosmological con- 
stant, or with a zero curvature. The horizontal line gives the 



Fig. 6. Di,s/Dos versus redshift for the I population in 
C10024-I-1654 and various cosmological models. The curves are 
strongly dependant upon a scaling factor which should be given by 
the accurate redshift of a large multiple arc. Here the horizontal lines 
give the value infererred from the observation by using an isother- 
mal sphere with core radius and the parameters of the best KKF 
model (full central line) corresponding to a change of 0.1 in the 
redshift population of the blue galaxies (top and bottom line). The 
(f^A = 0, n = 0) curve is the upper limit for any models with null 
cosmological constant. Therefore, they seems all ruled out by the 
observations. Low-JIa models are marginally compatible with the 
observations but would require very high redshift for the I galaxies. 
Best models are within the range 0.6 < Qa < 0.85. 



observed ratio for the last critical line of the / population in 
C10024-I-1654 {Ri = 60"), computed with the best param- 
eters of KKF (fT = ISOOfcm/.sec). We see immediately that 
many models are rejected. In particulax there is no model with 
a null cosmological constant compatible with the data. This 
seems to rule out any open universe and a {0.m = l)-universe 
as well. On the other hand, models with cosmological con- 
stant larger than 0.65 are acceptable. The upper limit for the 
cosmological constant is not provided by the lens, but some 
conservative upper bounds are given from the fim inferred by 
the CNOC dynamical analysis of rich clusters (Carlberg 1996) 
which give Sim = 0.16 ± 0.2. For a flat universe, this gives 
0.65 < 0,A < 0.85 . Furthermore, if we consider the robust 
upper limit of Kochanek (1996) coming from the statistic of 
gravitationally lensed QSOs, the possible range is apparently 
converging on an almost unique value: Qa ~ 0.65. 

The redshift range permitted for the distant / population 
is between 2.5 and 6.5, in agreement with the previous rela- 
tive scaling method. The redshift of the population at the first 
critical line observed at the beginning of the dip is still found 
to be near 0.9 ±0.1 as expected from the discussion in 3 and 
4.1. 

In this method, the determination of Ri is the critical pa- 
rameter. The uncertainties are dominated by the lens model, 
particularly by the error on the velocity dispersion. Assuming 
a redshift between a = 2 to 2 = 1 for the giant arc, the KKF's 
best model gives a between 1200 and 1400 km/sec, whereas 
Wallington et al. (1995) give 1100 to 1400 km/sec. This uncer- 
tainty will be suppressed as soon as the redshift of the triple 
arc will be well known (an observational priority). Both de- 
terminations are in good agreement with the galaxy velocity 
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dispersion directly measured by Dressier et al. (1260 km/sec; 
1985). In fact, it is remarkable that Bonnet ct al. have already 
shown that the shear profile in C10024+1654 is compatible with 
a isothermal sphere with a = ViOQkm/sec and Zsources = 0.9. 
If the velocity dispersion we use from the lens modelling is cor- 
rect, then the uncertainty only comes from the rather small 
rcdshift range wo found for the population at Rmin which can 
be converted to an equivalent small error on the velocity dis- 
persion. For Szmin ~ 0.1, 5(7 is as low as 50 km/sec. On figure 7 
we have also plotted the observed ratio Dls/Dos for a = 1350 
km/sec and a = 1250 km/sec. One can see that the constrains 
are still reasonnable, but it clearly emphasises that the depen- 
dancc on a may be an important shortcoming. 

Waiting for a precise redshift measurement of the large 
arc in C10024-I-1654 arc one have tentatively try to check the 
consistency of the results by searching the last critical line of 
the / population on another cluster for which the redshift of 
large multiple arc is known; A370. 




Fig. 7. CCD image of A370 showing the delimited areas where the 
counts have been performed. These are regions where crowding by 
bright galaxies is small and where no additional mass perturbation 
is detected. 



4-3. Scaling the last critical line in A370 

The redshift of the arc is 0.725 (Soucail et al. 1988) and 
the galaxy velocity dispersion is 1300 km/sec (Mollicr et al. 
1988). The cluster is at a redshift of 0.375, almost identi- 
cal to C10024-I-1654. Unfortunately, A370 is not as regular as 
C10024-I-1654. Its bimodal structure is well established from 
the lens modelling (Kneib et al. 1993) and from the ROSAT 
X-ray map (Fort and Mellier 1994). Furthermore, the arclet 
distribution shows that there is probably another mass con- 
densation to the cast (Kneib et al 1994). We restrict our study 
to two 60°-sectors along the axis passing through the center 
of the two cD galaxies which appear almost free of additional 
perturbations and are not crowded by bright foreground galax- 
ies. They are shown on figure 9. Since we only have very deep 
/ images with very good seeing, we concentrate on this color 
only. A detailed analysis will be given in paper II and we just 
give here a preliminary analysis that uses the symmetry of the 
bi-modal potential. 

The two clumps have a remarkable alignement and almost 
similar geometrical and dynamical parameters (Kneib ct a. 
1993). The depletion curves can be rcfered to the baryccntcr 
of the two clumps, though the bimodal shape of A370 makes 
the relation between the last critical line and the redshifts of 
sources less obvious than C10024-I-1654. Since we have selected 
two symmetric areas aligned with the main axis of A370, we 
expect the critical lines to be almost the same in the opposite 
North and South sector cones. The depletion curves are plot- 
ted on figure 10. They are very noisy because the number of 
detected galaxies is much lower than around C10024-fl654. In- 
deed they do not display a broad minimum so clearly, though 
their shape reveals a depleted area similar to C10024-I-1654. We 
find the possible radius of the extreme critical line along the 
main axis to be at 60±5 arcsec, surprisingly similar to the value 
obtained for C10024-I-1654. The vertical line on figure 10 indi- 
cates the position of the giant arc critical line Rb(z = 0.725) 
along the main axis. 

We applied the same method as in C10024-I-1654 to find the 
cosmologies compatible with the data. The two cD potentials 
are identical, with the same core radii rc and velocity disper- 
sions a. Hence, along the main axis, the potential well can be 
written: 

n^,y) = ^"'^J'"^ (Vl + (l-e)(x-:Do)2 + t,2 + 

^/l + {l-e){x + xoy+y^'^ (6) 

where e is the cUipticity of the two clumps (identical), and 
(xo, yo = 0) gives their position with respect to the barycenter. 

From '^{x,y) one can derive the position of the tangen- 
tial critical lines along the main axis. This gives a relation 
between the velocity dispersion, core radius and scaling coefH- 
cient Dls/Dos 

Dls_Qttg^ ( 1 ^ 

Dos rcc2 ^y'l-|-(l-e)(a;-a;o)2 



^\ + {\-e){x + xoY 
which can be solved numerically. 
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The term Qncr^ /rc(? can be obtained directly by using 
equation (7) at the position and the redshift of the giant arc. 
With a simple model assuming that e does not change with 
radius, equation (7) can be plotted in order to find acceptable 
solutions for the last critical radius observed in A370. 

If it is true, we found similar constraints on the redshift of 
sources and on the cosmological parameters as in C10024+1654. 
The most distant / galaxies corresponding to the last critical 
radius have 2.7 < z < 5 whatever the cosmological model. 
More precisely, if 0„i is larger than 0.16 from the CNOC works 
(resp. 0.3 from weak lensing in clusters) then the / galaxies are 
at a redshift lower than 6 (resp. lower than 4.7). 

Models with zero cosmological constant and Q.m < 0.15 are 
incompatible with the data. However, wo cannot rule out any 
values for Oa, even those close to 1. A similar remark has been 
incidently made by Kneib et al. (1993) when trying to explain 
the large ellipticity of the arclet A5 in A370 with the lens mod- 
elling. The uncertainty probably comes from the general shape 
of the cluster which is not as regular as our analytical descrip- 
tion. Fort & Mellier (1994) show the isoluminosity contours 
of the red cluster galaxies and the X-ray map. If light traces 
mass, then the discrepancy with respect to an regular ellipti- 
cal potential is clear. In this case, it is more appropriate to 
map the magnification from the local depletion N{x, y,< rn) 
everywhere in the cluster field. This preliminary study of A370 
confirms that the method of the last critical line gives reliable 
result only if it is used in conjonction with a self consistent 
modelling of the large arcs (paper II). 

5. Discussion 

The analysis of the magnification bias proposed by Broadhurst 
(1995) around rich lensing clusters can be used to infer the red- 
shift distribution of faint galaxies. First tentatives in C10024 
and A370 show that consistent informations were provided by 
these two clusters. The beginning of the dip at the first critical 
line corresponds to a redshift of 0.95±0.1 for the B and / galax- 
ies. The magnified / galaxies spread over a larger redshift range 
from ^ 0.9 to probably > 4 in both C10024 and A370. From the 
study of C10024-I-1654 we suggest that about 20% of the faint 
/ population is at redshifts above z = 4, and together with the 
knowledge of the lens model, namely the velocity dispersion, 
they give interesting information on the cosmology. It is sur- 
prising to find that both for the B and / selected galaxies the 
abondance of galaxies is not a uniformly decreasing function 
of redshift but rather corresponds to a bimodal distribution 
with few galaxies in the redshift bin 1.3 ~ 2. The shape of the 
distribution is remarkably similar to the photometric redshift 
distribution of the HDF population predicted by Stephen ot al. 
1996. Our sample cannot reproduce the shape below z = 0.4 
for obvious reasons, but we clearly have a peak below z = 1.1 
and a secondary broad peak between z = 1.5 to z > 3. The 
apparent bump in the middle of the depletion curve in I could 
be real bias magnification signature but it could also be related 
to the presence of an additionnal condensation of galaxies. 

The use of the velocity dispersion obtained from the best 
lens modelling of C10024-I-1654 gives the position of the in- 
ner critical lines with a redshift consistent with our findings 
from the spectro-photomctric observation of blue arc(let)s. The 
spectrum of the C10024-I-1654 arc is also a preliminary confir- 
mation of the redshift prediction. For the most outer critical 
line the significant fraction of the red galaxies which are not 
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Fig. 8. Depletion curve similar as figure 2 and 3 obtained from 
tlie deep I counts in A370. Tiie dotted line show the count for the 
notliern region only and the dashed line for the southern cone. The 
full line is the addition of the two regions. The vertical line show the 
position of the giant arc which provide and accurate image of the 
blue critical line. The depletion is clearly visible and shows the same 
characteristic as the I curve for C10014-I-1654, though the signal to 
noise ratio is smaller. In particular, the last critical line is 60 ± 3 
arcsec. 



detected in the S-band suggests that the Lya break lies be- 
yond 5000A, which also favors a redshift larger than 3. These 
galaxies might be the faint end of population detected with 
the Keck telescope by Steidel et al. (1996). Prom table 1 there 

number density is about 4 ± 3 / arcmin^ . 

For the cosmological test with the last critical line, a dif- 
ficulty when using the magnification bias from galaxy num- 
ber counts is the clustering of the background galaxies that 
may adulterated the curve N{r). Highly clustered background 
sources may bias the depletion curves toward a well-defined 
peak at some radius. It seems this major difficulty has been 
mostly avoided by going to fainter background galaxies. An 
additionnal test of the surface brightness to select the most 
distant ones with less clustering (see section 3) is too difficult 
for objects close to the noise. The fact that the two curves 
obtained here fit well those expected from the lens modelling 
seems an a-posteriori confirmation that the method is efficient. 
Besides, a large condensation of blue galaxies just behind the 
cluster center is not visible as a bump in the galaxies count in 
the cluster field with respect to our blank field or to the HDF 
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Fig. 9. Ratio Rj/Rb versus redsliift of tlie / sources in A370 for vEir- 
ious flat universes with cosmological constant. The scaUng is given 
by the redshift of the giant arc (0.724), so we do not need to plot 
curves for various redshifts as for C10024 (see fig. 5). The horizon- 
tal line gives the same quantity as inferred from the observations of 
Rl = 60 and by using the simple model provided by equation (7). 
There is no model ruled out by these observations, but they show 
again that the redshifts of the / population observed on the last 
critical line should be larger in the range 3. — 4.5. 



Fig. 10. Same plot as figure 9, for (Ha = 0)-universes. The red- 
shift of the / population observed on the last critical line has 
a well-defined lower limit at z = 3., but the upper limit is not 
well constrained. However, redshift larger than 5 are unlikely since 
Qm > 0.15 from the mass estimates of clusters from the CNOC sur- 
vey. Prom the mass estimates obtained by strong and weaJc lensing 
analyses in cluster, Qm > 0.3 and we can reasonably conclude that 
z < 4.7. 



field. Figure 1 does not show any perturbation of the slope in 
the cluster field. 

The constraints we obtained for the cosmological param- 
eters are essentially based on the dctcrmiuation of the "last" 
critical lines in /. The possibility to materialize a critical line 
at very large redshift seems to provide a lower limit on the 
cosmogical constant since the deviation angle is larger than ex- 
pected for an open universe. Some effects may affect the results. 
One is the fact that wc arc observing clusters with large arcs 
which could be sources associated with a large condensation 
of mass also. This mass might act as an additional deflector 
for most distant I galaxies that would then be deflected twice 
when the light passes through the two consecutive lenses. This 
might be wrongly modelled as a single lens with a larger critical 
radius. With deeper exposure time with a VLT a test can be 
performed to detect such possible condensation of galaxies by 
using color index distribution of galaxies (Mellier et al. 1994). 
The frequent discovery of unexpected shear around foreground 
groups in the field of over-luminous quasars (Bonnet ct al. 1993, 
Fort et al. 1996) is an indication that mass condensations are 
not rare. Similarly the "elliptical crisis" that come from the too 
high occurence of four lensed quasars can be also explained if 
their field has an additional shear field. These additional defiec- 
tors arc indeed particular cases of additional mass sheets which 
increase the angular scale of the gravitational bench as a cos- 
mological constant would do. Such intervening matter could 
be difficult to observe since it does not modify the shear pat- 
tern (Schneider & Seitz 1995). A second possible bias is that 
the ellipticity of the cluster potential changes significantly be- 
tween 30' and 60" from the cluster center but this is not seen 
on the shape of the depleted area around C10024+1654. Al- 
though these two bigises cannot be systematic for all cluster 
with large arcs they have to be kept in mind for C10024-I-1654 
and A370 which display several bright blue arclets. With the 



lens model it is possible to infer their clustering properties in 
the source plan if we suppose they are all at the same redshift. 
This should be done in the future. 

In conclusion, going close to the noise to search for an 
extreme critical line can provide information on the redshift 
distributions of the faintest galaxies and a direct test on the 
cosmology. It seems from this preliminary investigation that a 
flat universe with a cosmological constant near 0.65 could be a 
possibility. This value is still in agreement with the Kochanek's 
(1995, 1996) upper value coming from the statistical analyses 
of lensed quasars, and from the tentative predictions of Mo 
and Fukugita (1996). It is also compatible with the parame- 
ter range inferred by van Waerbeke et al. (1996b) from the 
neo-V/Vmax statistics of quasars. As well, the result for a large 
abundance of high-redshift galaxies is compatible with the re- 
cent spectroscopic detection with the Keck telescope by Steidel 
et al (1996). An interesting consequence of the large radial dis- 
tance of the last / critical line is that red gravitational pairs 
candidates should be now searched for much larger radii than 
the large blue arcs. A few candidates are known around some 
clusters and their spectroscopy could give a conflrmation of all 
the results presented here. It is also clear that it is absolutely 
necessary to use our method on a large numbers of clusters 
and to use very good modelling of large gravitational arcs sys- 
tems before to expect strong constraints on the cosmological 
parameters. Such an observing programme is a typical one for 
10-meters telescopes. 
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